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Abstract. We provide a leading order semiclassical asymptotics of the 
energy of bound states for magnetic Neumann Schrodinger operators 
in two dimensional (exterior) domains with smooth boundaries. The 
asymptotics is valid all the way up to the bottom of the essential spec- 
trum. When the spectral parameter is varied near the value where bound 
states become allowed in the interior of the domain, we show that the 
energy has a boundary and a bulk component. The estimates rely on 
coherent states, in particular on the construction of 'boundary coherent 
states', and magnetic Lieb-Thirring estimates. 



1. Introduction 

Let Q' C M? be an open and bounded domain with regular boundary. We 
will consider both the case of interior domains = 17' and exterior domains 

= r2 \ n'. 

Consider the magnetic Schrodinger operator, or magnetic Laplacian, in 
It has been observed by many authors (see for example [l2]) that the presence 
of Neumann boundary conditions has an effect similar to that of a negative 
electric potential (also Robin boundary conditions have a similar effect, see 
[T3j[ll]). For the case of the present discussion, consider a constant magnetic 
field, more generally we will impose the hypothesis (|1.6|) below. One aspect 
of the analogy is that the Neumann boundary condition leads to a discrete 
spectrum below the lowest Landau level (in the case of an exterior domain, 
the (magnetic) Laplacian may have essential spectrum). It is that spectrum 
which we will discuss in the present paper. 

One motivation for our work is the analysis of type II superconductivity 
for applied magnetic fields close to the second critical field (see [T71 Prob- 
lem 2.2.8, p. 491] for a discussion of this point). In the setting of supercon- 
ductivity, one encounters the asymptotic regime of a large magnetic field, but 
this is equivalent (through a simple change of parameter) to the semiclassical 
regime considered in this paper. 

Many different investigations of the discrete spectrum close to the lowest 
Landau level have appeared (see [TOl [191 EO] and references therein). Of 
particular importance for our present investigation is the recent paper [9]. 
In that paper the counting function is studied up to an energy strictly below 
the lowest Landau level. Inspired by work on the magnetic problems without 
boundaries — but with electric potential — see [16], we shift focus from the 
counting function to the energy, i.e. the sum of the eigenvalues. This allows 
us to obtain 'semiclassical' results all the way up to the lowest Landau level 
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(Theorem 11.11 below). Furthermore, by varying slightly the energy up to 
which we sum the eigenvalues, we can demonstrate how the bulk or interior 
bound states start to contribute to the leading order asymptotics for the 
energy (Theorem 11.41 below). 

Let us now move to the precise statement of our results. We consider the 
(Neumann) Schrodinger operator with magnetic field : 

= -(/iV-iA)2 = (-z/iV- A)2, (1.1) 

whose domain is, 

D {Ph,n) = {u £ L\n) : {hV - iAy u £ L^Q), j = 1,2, (1.2) 

u • (W - iA)u = on dn}. 

Here /i > is a small parameter (the so called semi-classical parameter), u 
is the unit outward normal vector of the boundary dQ and A G C'^{^l;M?) is 
a given vector field — the vector potential. The magnetic field is a function 
and is given by 

B = curlA = d^,A2-d^,Ai. (1.3) 
With this magnetic field we associate the quantities 

b=mi_B{x), b':=mfB{x). (1.4) 

Assuming that 6 > 0, we know that the following leading order asymptotic 
expansion holds for the bottom of the spectrum of Ph n (see for example 

m) 

inf SpecP/i^n = hmm{b,Qob') + o(/i). (1.5) 

Here Go S]0, 1[ is a universal constant (the definition will be recalled in 
(I2l3]l below). 

We shall assume that the magnetic field is bounded, positive and satisfies 

b > @ob' > 0. (1.6) 

Notice that the hypothesis p.6p is satisfied in the case of a constant magnetic 
field. Under the hypothesis p.6p eigenvalues strictly below bh are associated 
with eigenfunctions localized near the boundary. One of the objectives of 
the present paper is to prove that — for the energy — this remains true all the 
way up to the value bh. 

In order to state our main results, we need to recall some facts concerning 
the harmonic oscillator on the semi-axis M+. For ^ G M, we denote by /ii(^) 
the lowest eigenvalue of the operator 

-df + {t-Cf in L\R+) 

with standard Neumann boundary condition at the origin. It is well known 
(see OH]) that the function i— > ^i(^) is smooth, 

/ii(C) < 1, for all C G M+, MO > 1, for all ^ G M_, (1.7) 
and the integral 

/ (^i(e)-i)dc = - iMo-M-d^ 

Jo JR 
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is negative and finite. Here we introduced the notation more generally, 
we will use the following positive functions 

j X, X > 0, Jo, X > 0, 

X + = < X _ = < 

[0, X < 0, |-x, X < 0. 

Our result is the following. 

Theorem 1.1. Suppose is either an exterior or an interior domain. Sup- 
pose B satisfies p.6p . Then the spectrum of Ph,n below bh is discrete, 

(^{Ph,n)n] - (x,bh[= {ei{h),e2{h),...}, 

and the sequence {ej{h) — bh}j>i is summahle, 

XI [^i(^) ~ ^^]- = {{Ph,Q. - bh)l]_^^hh[{Ph,n)) is finite. 
j 

Moreover, the following asymptotic formula holds. 



^- [ r B{xf/'\-^-f^,{d dCds{x). (1.8) 



27r 

Here ds{x) denotes integration with respect to arc-length along the boundary 

dn. 

Remark 1.2. Assuming that the magnetic field B = curl A is constant, 
B{x) = b, the asymptotic formula of Theorem 11.11 reads 

Ihn [e,ih) - bh]_^ = (1 - ^i(C)) d^ . 

Remark 1.3. In [9], an asymptotic formula is obtained for the number 
N{Xh) of eigenvalues of Ph^n below A/i, for a given X < b. The precise result 
is the following : 

lim h^/^N{Xh) = ^ [ B{x)^'^ dids{x) . (1.9) 

Integration of (|1.9p yields the following formula for the energy. 



\iu.h-^/^y^e^{h)-Xhl 

n — >(j 



1 

2^ 



an Jo 



S(x)3/2 



B{x) 



d^ds{x). (1.10) 



However, the proof we give to Theorem 11.11 gives equally (11.101) (we only give 
the details for the harder case A = 6), and hence, by differentiating (11.101) . 
(more precisely, the 'differentiation' needed is the technique used to go from 
energies to densities in semiclassical problems, see [6] for details) we provide 
an alternative proof of (II. 9p . 
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In the next theorem, we restrict to the case of interior domains, i.e. 
bounded In this specific case, the operator Ph^n has compact resolvent 
and hence its spectrum is purely discrete. Let us denote by {ej{h)} the 
increasing sequence of eigenvalues of Ph,n (counted with multiplicity). 

Theorem 1.4. Suppose S7 is bounded, smooth and B is constant in VL. Given 
a G M, the following limit holds, 



Remark 1.5. The first term on the r.h.s. of p.lljl corresponds to boundary 
states (i.e. eigenfunctions localized near the boundary) and the second one 
corresponds to bulk states. Theorem 11.41 sharpens our understanding of the 
transition from boundary states to bulk states. 

Remark 1.6. Since inf SpeCj^gg Ph,Q, = bh in the case of exterior domains, we 
see that the left hand side of (|1.11|) diverges for any a > 0. That forces us 
to restrict to bounded Q.. Our methods would also apply to non-constant B, 
but the order to which the bulk term appears depends on the local behavior 
of B near the set {x : B{x) = b}. For simplicity of exposition we therefore 
restrict to the case of constant field. 

The approach we follow is inspired by that of Lieb, Solovej and Yngvason 
[16| . but extended with the construction of coherent states for the half-plane 
operator. In order to control the errors resulting from the approximations, 
we make use of a Lieb-Thirring inequality for magnetic operators (see [5l [E] 
and references therein), and a rough estimate of the energy of bound states 
for the case of a cylindric domain. 

The paper is organized as follows. We collect in Section [2] some prelimi- 
naries. In Section EJ we determine a rough bound for the energy of bound 
states of the operator p.ip in the case when the domain is a cylinder. 
Section m is devoted to the construction of coherent states. In Section we 
prove Theorem II. 1[ Finally, in Section [6] we give the additional details to 
achieve the proof of Theorem II. 4[ 



2.1. Lieb-Thirring inequality. Let B £ C^{R'^;R) n L°°(M?) be a mag- 
netic field such that B(x) > for all x £ M^. The vector field defined 
by 



We have the following Lieb-Thirring estimate for the negative eigenvalues of 
H^2 + V (see [5] and references therein). 





2. Preliminaries 




provides a magnetic potential for B = curl^. 
Consider the Schrodinger (Pauli) operator 

H^^ = -(y -iAf - B mL^iR^). 
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Theorem 2.1. There exists a universal constant C > such that the fol- 
lowing estimate is valid for the sum of the negative eigenvalues {ej}j>i of 
the operator H = H^2 + V , 

Y,\e,\<c(\\Bh^^u^) [ [V].dx+ [ [Vtdx 

J \ JR2 Jr2 

2.2. Variational principles. Let H he a self-adjoint operator in L^(M^) 
(of domain D{H)) such that 

inf C7ess(^^) > 

H1]_^Q[{H) is trace class. 

We shall need the following two simple variational principles concerning 
the operator H, which are frequently used in [TBI [5]. 

Lemma 2.2. Let ^ be a bounded operator such that < 7 < 1 (in the sense 
of quadratic forms) and the operator is trace class. Then it holds that, 

tr(^l]-oo,o[(^)) <tr(//7)- 

Lemma 2.3. Assume that the operator H satisfies the hypothesis (H). Then 
it holds that, 

N 

tr (m]_^,o[(^)) = inf , Hf,) , 

where the infimum is taken over all orthonormal families {/i, /2, . . . , /at} C 
D{H) andN> 1. 

2.3. Boundary coordinates. The closed quadratic form associated with 
the operator (jl.lj) is 



/ \{hV - iA)u\'^ dx , (2.1) 
Jn 



qh,n[u) 

with the form domain 

H^Ai^) = {u£ L^i^) : {hV - iA)u G L^{^)} . (2.2) 
The following magnetic potential generates a constant unit magnetic field, 

Ao(a;i,X2) = (-2:2,0), V (xi, X2) € M x M . (2.3) 
The quadratic form, 

H^^j^^^iR X R+) 3 u ^ \{hV - ibAo)u\'^ dx , (2.4) 

with 6 > 0, will serve as a model form. Actually, when the function u G 
-^/ia(^) supported near the boundary, the form (|2.4|) turns out to ap- 
proximate ()2.1I) . In order to make this precise, we introduce a convenient 
coordinate transformation valid in a sufficiently thin tubular neighborhood 
of the boundarjQ. For more details on these coordinates, see for instance [T2l 
Appendix A]. 

^to ■■ ^{to) 3 X ^ {s{x),t{x)) G X ]0, tol, (2.5) 



^Here we assume for simplicity that dO, is connected. In general, dO, has finite connected 
components and therefore we should work on each component independently. 
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where for Iq > 0, 0,{to) is the tubular neighborhood of d^l: 

n{to)={x£n : dist{x,dn) <to}. 

Let us mention that t{x) = dist(x, dCl) measures the distance to the bound- 
ary and s{x) measures the curvihnear distance in dQ. We shall use the 
usual identification between the circle and the interval [0, \d^l\ [. The 

Jacobian of the transformation is equal to 

1 - tk{s) , 

where k denotes the curvature of 557. 

Since to will be fixed once and for all, we will sometimes omit it from the 
notation and simply write $ instead of ■ 

Using the coordinate transformation we associate to any function 
u e L'^i^), a function u defined in [0, \d^\[ x [0, to] by, 

uis,t) = u{<^^^\s,t)). (2.6) 

Furthermore, the function u extends naturally to a |9r2|-periodic function in 
sGR. 

We get then the following change of variable formulae. 
Proposition 2.4. Let u G H]^{i}{tQ)). We write u{s,t) = u{^tois,t)), 
Ai=Aio$i^, A2 = A2 0$i„. 

Then we have : 

r r\d^\ /■*() r 

/ \{V - iA)u\^ dx = / / \{ds-iAi)u\ 
Jn(tn) Jo Jo '- 



1 2 



-2|/Q -A \~\2 



+ (l-t/c(s))-^|(9t-zA2)Sr {l-tk{s))dsdt, (2.7) 
and 

\u{x)\^dx= / \u{s,t)\^{l-tk{s))dsdt. (2.8) 



/Q{io) 

We will use the symbol C/# for the operator that maps u to u. We shall 
frequently make use of the next standard lemma, taken from [9], Lemma 3.5]. 

Lemma 2.5. There exists a constant C > and for all Si G [0, |(9r2| [, 
52 €]5'i,|9r2| [, there exists a function (p S Cq{[Si, S2] x [0,to];M) such that, 
for all 

Se[Si,S2], Te]OM, ee[CT,Cto], 

and for all u G Hj^p^{n) satisfying 

suppuC [81,82] X [0, T], 
one has the following estimate, 

qh,n{u) - I I (W - iBAo)e'^/^u\^ dsdt 

JrxR+ 

< [ (e\ (W - iBAo)e''^u\'^ + Ce"^ ((5^ + T^f + h^) dsdt. 

Here, 8 = 82 — 81, B = B{8,Q), the function u is associated to u by ^ 
and extended by onRx \ supp u. 
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2.4. A family of one-dimensional differential operators. Let us recall 
the main results obtained in [U [12] concerning the family of harmonic os- 
cillators with Neumann boundary condition. Given E R, we define the 
quadratic form, 

B\R+)Bu^q[C]iu)= [ \u'{t)\'^ + \{t-Ou{t)\'^dt, (2.9) 

where, for a positive integer G N and a given interval / CM, the space 
B^{I) is defined by : 

B''{I) = {ueH^{I)] fu{t) e L'^il), yj = l,...,k}. (2.10) 



Since the quadratic form p.9p is closed and symmetric it defines a unique 
self-adjoint operator C[C]- This operator has domain, 

D{£[C]) = {u e B\R+y, n'(0) = 0}, 

and is the realization of the differential operator, 

m = -df + it-o', (2.11) 

on the given domain. We denote by {/Uj(0}j=i the increasing sequence of 
eigenvalues of which are all simple. By the min-max principle, we have, 

,^iO= inf /l^. (2.12) 

It follows from analytic perturbation theory (see [15]) that the functions 

are analytic. 

As recalled in the introduction, /ii(0) = 1 and |/Ui(^) — 1| decays like 
exp(— ^^) as ^ — > -|-oo (see [3]) thus yielding that 

/ (^i(e) -l)d^ = - [/xi(e) - 1]- is finite. 
Jo Jr 

We define the constant : 

@o = mlfii{0. (2.13) 



Let us recall an important consequence of standard Sturm-Liouville theory 
(c.f.. [91 Lemma 2.1]). 

Lemma 2.6. The second eigenvalue satisfies, 

infMe)>l. 

Notice that part of this conclusion is a consequence of the analysis of 
Dauge-Helffer [4], who show that the infimum of ^2(0 is attained for a 
unique .^2 G R- 
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3. Rough energy bound for the cylinder 

Let us consider the operator (jl.ip in the particular case of a cyhndric 
domain 

n = [0,S]x]0,h^/'^T[. 

Functions in the domain of Ph^n satisfy Neumann condition at t = 0, periodic 
conditions at s G {0, 5} and Dirichlet condition at t = h^^'^T. We assume 
in addition that the magnetic field is constant B{x) = 6, 6 > 0, and that 
the magnetic potential A = Aq is the one given in (12. 3p . In this particular 
case, the operator has compact resolvent, hence the spectrum consists of an 
increasing sequence of eigenvalues (ej)j>i converging to +oo. In particular, 
given A > 0, the energy 

£iX,b,S,T) = Y,[hb{l + X)-ej]^ (3.1) 
j 

is finite. Our aim in this section is to provide a rough estimate of this energy. 

Lemma 3.1. There exist positive constants Tq and Aq such that, for all 
S > 0, b > 0, T > VbTo and A s]0, Aq], we have, 

£{X, b, S, T) < (1 + X)hb + 1 

Proof. By separation of variables and a scaling we may decompose Ph^Q as 
a direct sum: 

^hb(-^ + i2nnh'/H~'/'S~' + tf) in ^ L\]0,T/Vb[) , 

with Neumann boundary condition at the origin and Dirichlet condition at 
t = T l\fb. Therefore, we may express the energy (|3.1|) in the form, 

f(A,6,5,T) = /i6^ [l + A-;Uj(27rn/i^/2^-^/25-i.2-/^) (3 2) 

Here, for a given ^ G R, we denote by Hj{^;T) the increasing sequence of 
eigenvalues for the operator 

with Neumann condition at the origin, and Dirichlet condition at t = T. No- 
tice that for a fixed T, the min-max principle gives immediately, /ii(^;T) — > 
+00 as ^ — > +00. Hence, the sum on the right hand side of (13. 2|) is finite. It 
follows also from the min-max principle that //2(C;'^) ^ /^2(0: where //2(0 
is the second eigenvalue of the operator (12. lip . Now, Lemma [2?6l gives the 
existence of a sufficiently small Aq such that /i2(0 > 1 + Aq for all ^ G M. 
Thus, taking A g]0, Ao[, 

£{X,b,S,T) =hb^^l + X- fii{2TTnh^/^b~^/^S-^;T/Vb) . (3.3) 

By (|3.3|] . if one can localize the set, 

{eeM : M^;T/Vb) < 1 + A}, 
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then one gets immediately an estimate of the energy £{X,b, S,T). Notice 
that for t < T and |^| > 2T, it holds that {t - > , hence by the min- 
max principle, /ii(^;T) > . Therefore, choosing Tq > 1 + Aq, it holds for 
T > 7^3 and A G]0, Ao[, 

fii{^;T) < 1 + \ =^ lei <2r. 

From the above localization, the estimate of Lemma 13.11 becomes a conse- 
quence of (|3.3I) . □ 

4. Generalized eigenprojectors 

4.1. Eigenprojectors in M^. Let Aq be the magnetic potential given in 
(|2.3p . b > and Ph,b,R^ the self-adjoint operator 

Ph,h,RxR+ = -{hV - ibAof in L\R^) . 

We recall in this section the well-known eigenprojectors for the operator 
P/jb]K2. These are projectors Ilj{h,b), j > 1, on the Landau levels that 
satisfy in particular, 

oo 

Ph^,^^2Uf{h, b) = (2j - m'-ih, b) , nj"!^, b) = Idi2(IR2;C) • 

i=i 

The integral kernels (denoted again by n^(/i, b)) are given explicitly (see [IH 
(3.13)]), 



n^(/i,6)(2;,y) = TT^exp 



ib{yiy2 - X1X2) 



2TTh ^ V 2/1 

f ib{xiy2 - X2yi) b\x-y\'^\ f b\x - y\ 
X exp ( — — ) Lj-i 



2h 4h J ' \ 2h 

where Lj are Laguerre polynomials normalized so that Lj{0) = 1. We will 
need the fact that, 

Ilf{x,x) = ^, VxeM^ (4.1) 

4.2. Eigenprojectors in the half-space. We construct in this section pro- 
jections on the (generalized) eigenfunctions for the operator (II. 1|) in the case 
= R X M+. 

Let Aq be the magnetic potential given in (|2.3p . 5 > and Ph,b,RxR+ the 
Neumann realization of the operator 

Ph,b,RxR+ = -{hV- ibAof in L^{R x R+) . 

Let us denote by (uj(-;0)jli orthonormal family of real-valued eigen- 
functions of the operator C[^] from (|2.1ip . i.e. 



-n;'(t;0 + {t- e)S (t;0 = /^i(6%(i;0 , in M+ , 

/ Uj{t;if dt = l. 
Jr+ 

Let us define a bounded function M x IR_|_ 3 (s,t) 1— > Vj{s,t;(,) by : 

Vj (s, t ; = exp (-ie s) uj {t; ■ (4.2) 
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We introduce a family of projectors nj(^) on the functions Vj , 



C^(]R X M+) 9 / Vj{xi,X2;0vj{yi,y2;C) <fiyi,y2)dyidy2. (4.3) 

The projectors 11^ are kernel operators. For j £ N, the integral kernels are 
defined by: 



Kj{0{ixi,X2),{yi,y2)) = Vj{xi,X2;^.) x Vj{yi,y2;0 ■ (4-4) 
Now, one easily verifies the following properties : 



oo „ 

,=1 



id^ = 27rIdi2(KxM+) , (4.5) 
and 

Pi,i,RxR+nj(e) =^,(Onj(0- (4.6) 

By means of a dilation, we get a family of eigenprojectors for the operator 
Ph,b,RxM.+ - Let US introduce the unitary operator, 

Uh,b : i^(M X IR+) 9 ^ Uh,bV G L'^(R x M+) , (4.7) 

such that, for all x = {xi,X2) G M x 

iUh,b^) (x) = y'li/h ip{y'l)/h x) . 

Notice that 

Uh,b-^h,b,RxR+Uh,b = /i&-Pl,l,RxR+ • (4.8) 

Then we introduce the family of projectors, 

Uj{h,b;0 = Uh,bIl,{OU^i. (4.9) 

Again, the projectors Ilj{h, 6; ^) are kernel operators. For j £ N, the integral 
kernels are given via the kernels ()4.4p : 



K,{h,b;0{x,y) = j^KjiO {Vb/hx,^/b/hy) 



= ^e-'V^^^^^-y^^u,{^/b/hx2-,0ujiVb/hy2-,0 , (4.10) 

for all x,y G M x IR_|.. Now, the following properties are directly inferred 
from KB and KEh : 



oo „ 

/ n,(/i,6;e)a!C = 27rIdi2(R^K+), (4-11) 

,=1 



and 

Ph,bMxR+'^3{h,b;0 = hbnj{^)Uj{h,b;^). (4.12) 
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5. Proof of Theorem 11.11 

5.1. Existence of discrete spectrum. Recall the operator P/^q intro- 
duced in (II. ip . Tlie following inequality holds (see [1]) for compactly sup- 
ported functions 

[ \{hV - iA)^\^ dx > h [ B{x)\(p\'^ dx , V (/? G C^(J^) . 
Jn Jn 

Using then a 'magnetic' version of Persson's Lemma (see OdH]), we get that 

inf Spec^ss Ph,n > bh , 

hence proving the first statement of Theorem 11.11 Let {ej{h)} be the se- 
quence of eigenvalues (counted with multiplicity) corresponding to (7{Ph,n)f^] — 
oo, bh[ , and define H = Ph,n — bh. It suffices now to show that 

^(/,,i7/,)>-oo, 

3 

where the functions fj are normalized eigenfunctions associated with the 
eigenvalues Cj. To that end, we introduce a partition of unity of M, 

V'i+'02 = 1> supp^i C] - OO, 1[, suppV'2 C [^,oo[, (5.1) 

and set for /c = 1,2, Xk{x) = ipk{t{x)), x G M^, where t{x) is the signed 
distance to 9^1, 

t{x) = dist{x,dQ) ii X £ , t{x) = —dist{x,dQ) otherwise. 
By the IMS formula, we write, 

2 

(/, , Hfj) = J2 {{Xkfj , Hxkfj) - \VxkWf) , (5.2) 

k=l 

where || • || denotes the norm in Q. Using the bound B > b together with 
the fact that ipf + tp2 = we get a further decomposition of (|5.2p . 

2 2 

if J , Hf,) = Y,{fj , Xk{H - V)xkfj) > Y^ifj , Xk{Ho - V)xkfj) , (5.3) 

k=l k=l 

where Hq = — (/iV - iA)^ - hB and 

^ = /i'(|Vxi|' + |Vx2|') . (5.4) 

Pick an arbitrary positive integer N < Card({ej}j). Let us define the trial 
density matrix L^{R'^) B f ^ -f2f ^ L^iU"^), 

N 

72/ = X2^(X2/, fj)fj, 
J=l 

which verifies the conditions < 72 < 1 (in the sense of quadratic forms) and 
{Hb — V^)72 is trace class (actually this is a finite-rank operator). Moreover, 
using Lemma [2711 we know that the operator {H^2 — V^)l]-oo,o[(-f^R2 — V) is 
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trace-class, where H^2 = —{hV — iA)'^ — hB is now the corresponding Pauh 
operator acting in L^(]R^). Therefore, we deduce from Lemma [2^ 

N 

Y,{fj , X2{Hb - V)x2fj) = tr[iHB - F)72] 

> tr {{H^2 - y)l]_oo,o[(^R2 - V)) , (5.5) 

and we notice that this bound is uniform in N. 

To take care of the boundary contribution, we set to = mt{Q n suppxi) 
and define the trial density matrix 

N 

71 : L^ioj) 9 / ^ 71/ = XI ^^(xi/ , f,) fj G L^{oj) ■ 

i=i 

Thereby, we also get, 

if J , xi{Hb - V)xifj) = tT[{HB - Vhi] 

>tr{{H^-V)l]_^^o[{H^-V)) , 

where H^^ is the restriction of Hb on uj with Dirichlet condition on 
U n (9(suppxi)- One should notice that, since u! is bounded, H^^ — V has 
compact resolvent and hence {Hi^ — l^)l]-oo,o[(-f^aj — V) is evidently trace- 
class. Coming back to (|5.2|) and (|5.3p . we deduce that. 



- Y.[e,ih) - bh]^ = Y,{fj , Hfj) > tr {{H^. - y)l]„^,o[(^R2 " V)) 
j j 

+ tr{{H^-V)l]_^^o[iH^-V)) , 
proving thus that the sequence ej{h) — bh is summable. 

5.2. Lower bound. Let {/i, /2, • • • , In} now be any orthonormal set in 
D{H). We will give a uniform lower bound to 

N 

i=i 

Using Lemma [231 this will imply a lower bound to tr (^H1^_^ q[{H)Y 
Step 1. Localization to the boundary. 

Let r(/i) g]0, 1[ be a small number to be chosen later. Using the partition of 
unity in (15. Ih . we put 



^i,/.(^) = V'l (^) , V'2,h(x)=V'2(^), vxea 

By the IMS formula, we write, 

2 

(/,• , Hfj) = J2{fj > - Vh)^k,hfj) , (5.6) 

fc=i 

where 

Vh = (|V^i,fc|2 + \Vi^2,h?) . (5.7) 
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Notice that the term corresponding to A; = 2 in (|5.6p corresponds to the 
interior term. We will prove that it is a lower order error term. Similarly to 
(|5.5p . one can show that, 

and we remind the reader that H^2 = — (/iV — iA)'^ — hB{x) is now the cor- 
responding self-adjoint operator acting in L^(R^). Using the Lieb-Thirring 
inequality of Lemma 12.11 we get 

N 



> -c 



[ {\\B\\L^[-h~Wh]- + i-h-Wht) dx] . 



With our potential Vh from ()5.7p . the integral on the right side above becomes 
of the order of + 7^^- Thus, we get 



N N 

h / h 

rjh) ^ V{hy^ 

Later, we shall choose t(/i) in such a manner that the first term on the right 

hand side above is the dominant term. 

Step 2. Boundary analysis. 

Here we consider the boundary term, 

N 

Y,{f,,^i,h{H-Vh)^i,hf,). 
i=i 

To that end, let x ^ be positive, smooth, supported in ]0,1[ and 

satisfying 

Jx\s)ds = l. (5.9) 
Using the boundary coordinates {s,t) introduced in (12. Sp . we put 

Y^r(/i) V J 
and we notice that, for all cr G ]R\] — t(/i), \dQ\ [, 

Xh{x]cr)=0 VxGil(to). 
Using again an IMS type decomposition, we write, 

N 

N 

= / yZ'^fi ^ '^'^^hXh{x;a){H -Wh)ipi,hXh{x]cr)) da , (5.10) 



Cz7rv(l + Z7Mf ) • (5-8) 
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where 

Wh = Vh + h^ f \VXh{x; da . (5.11) 
Let us denote by (<I>t(, is the coordinate change (12. 5p valid near the boundary) 
Uj,h{(y) = Uj^h{x;a) = il)i^h{x)xh{x;a)fj{x), B„ = 5(^>j~^(cr, 0)) . (5.12) 

Using Lemma 12.51 we get a constant C > and a scalar function (p = cj)^ 
such that, for all CT{h) < e <C 1, we have, 

AT 

if J , ^i,hXh{<T){H - Wh)^IJi,hXh{T)) >{l-e) 

N 

X Y,{^"^'''%,h{<y) , {Hh,B^,RxR+ - m)e''^-^%A'^))LHR^R+) , (5.13) 
i=i 

where 

Hh,B^,RxR+ = Ph,B„,RxR+ — bh = — (W - iB^Ao)'^ - bh , 
Aq is the potential introduced in (|2.3p . 

Wh = Wh + Ce~' (t(/i)4 + h^) , (5.14) 

and to a function v{x), we associate the function v{s,t) by means of (|2.6|) . 
The functions Uj^h{a) are naturally extended by in M x \ suppnj^/i(a). 

Later, we shall make a suitable choice of the parameter e. We use part of 
the kinetic energy i^h,B^,RxR+ to control the error resulting from W/j. Let 
5 G]0, 1[ be another parameter to be specified later; 5 will be chosen as a 
function of h. We decompose the previous sum in two , 

N 

Y,{e^'*'^'%,h{<y) , (^h,B.,RxM+ - m)e''*'^'%,h{<y))L^iRxR+) (5.15) 

N 

= {l-5) Y,{e''^^l%,n{o) , i?h,B.,MxR+e*'^''/'^S,, ^(a))z.2(KxR+) 

N 

+ 6Y{e'^^/%Ma),{Hh,B.,RxR+-S-'Wh)e'^^^%A^))LHRxR+)- 
i=i 

Let us estimate the last term on the right side. We take T = T{h) » 1, to 
be specified later, and we make the following choice of r(/i), 

r(/i) = h^/^T. 

Notice now that 

P"'Wft||ioc(iRxM+) < C6-^ {e-'h^T^ + /iT-2) < XB^h , (5.16) 
where we have defined 

A = Cb-^6-^ [e-^hT^ + T'"^) . (5.17) 
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Furthermore, with S = T{h) and the above notation, we define the operator 
7 on L2([0,5]x]0,/iV2r[), 

TV 

The next estimate is a simple apphcation of Proposition 12.41 

0<(/,7/) (5.18) 

<{l + CT{h)) [ Xh{x;afi;lf,{x)\{fo^~^) {x)\' dx < C'rih)-^ 



L2 • 

Thus, (7/^(fe)-i 7 is a density matrix. The operator 7 is constructed so that 
we may write the term we wish to estimate in the following form, 

N 

i=i 

= tr [{Hh,B,^RxR+ - S"^^^h)l] ■ 

Using the bound (|5.16p on the potential Wh and the variational principle of 
Lemma [221 we get, 

N 

>-C'T{hy^8{\B^,S,T). (5.19) 

Here we remind the reader that S = T{h) = h^^^T, the energy £{X, B^j, S, T) 
is introduced in (|3.ip . and C is a positive constant depending only on J7. 

In order to use the estimate given in Lemma ISTTl we make the following 
choice of e and 5, 

T = , e = hT' = h^^ , 6 = r-3/2 = /,3/i6 ^ 

so that the following conditions hold : 

h^^^T < e < 1 , A « 1 . 

Notice that this choice implies 

r(/i) = (5.20) 

Using now Lemma ISTTl we conclude that (recall that S = T{h) in our case 
which compensates the factor from (I5.18P ). 

N 

> -Ch^/^TBl/^ , 

and consequently, after integration w.r.t. cr g] — t(/i), \dCl\ [, 
N 

> -C6Th^/^ = -Ch^/^^ , (5.21) 
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thus obtaining an error of the order o(/i^/^). We point also out that with our 
choice, r(/i) = h^/^ , so that the error in (15.8(1 becomes small, of the order of 

Step 3. Estimating the leading order term. 

Now we estimate the first term on the right hand side of (|5.15p . We use the 
projectors lij constructed in Section [H In view of ()4.11|) and ()4.12|) . we have 
the following splitting, 



BJi 
2tt 



p=i 

oo 



,B„,MxI[ 



Since i?o- > b, Lemma EJ] gives that for p > 2, fJ,p{C) — > 0. Hence, we 
get for any function / in the domain of -fffc,Ba,lRxR+ > 



{f,Ui{h,B^;()f)d^. 



Therefore, defining 

N 

\ / L/ fMxM.-|-] 



we get that 



(5.22) 



N 



B„h 



> 



2ir 



W(6 



L2(RxR+) 



57v(a;Oo^^ (5.23) 



We estimate now the term on the right hand side above. We start first by 
estimating 5Ar(cr;^). Recalling the definition of Uj^h, and using again the 
coordinate transformation valid near the boundary, we get (using the choice 



e^^^/%,.(.),n,(/.,5.;e)e*/%,.(.)^^^^^^^^^ 

<{l + Ch'/''){fj,r{h,B^;a,Ofj)mn) ■ (5-24) 
Here V{h, B(^;cr,^) is a positive operator, which is given by, 

V{h, B^;a, = i'i,hXh(.<^) Vt^""^"'"" ni(/i, B.'i) e"'^^'^ U^A,hXh{^) , 

and the transformation [/$ is associated to the coordinate change intro- 
duced in (12. 5p . Since {/i, /2, . . . , /at} is an orthonormal family in L'^{Q), we 
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deduce that 

N 

Y,{fj,r{h,B^;a,Ofj)mn) 

3=1 

< tr {V{h,B^;a,0) 



B. 



^Xh{s-,a)\^\Mt)\'\Mh~'/'Bl/H;0\'a-tk{s))dsdt, 



where 



\Xhis-,a)\'\Mt)\'\Mh-'/'Bl/H-0\'dsdt 



X 



s — a 



r{h) 



ds [ \ui{h-^/^Bl/h;0\'^dt 



< r{h)-^ 

= h^/^B~^'\ 
and similarly, 

|xh(s;a)PlV'u(t)lVi(/i"'/'^y'i;OI't^(^)dsdt 

= 0{h^i^T{h)) = o{h'^/^) . 

Coming back to (|5.22l) and (I5.24p . we get that, 

0<SN{c7;0<^ + Oih~'/^). 

Implementing the above estimate in (I5.23p . we deduce the following lower 
bound, 

N 



h^l^Bl' 



> 



2tt 



/"i(e) 



B, 



di - 0{h^'^) 



and consequently, upon integrating with respect to a (recall that u.j^h{a) is 
supported in {0 < s - < t(/i) : s G [0, [ }), 



Hani ^1/2^3/2 



> 



2tt 



B^ 



dida - 0{h^/^) . (5.25) 



Recalling the definition of B^ = i?($^^^((T, 0)), the integral on the r.h.s. 
above is actually nothing but 



27r 



B{xfl^ 



m(0 



B{x) 



d^\ dx . 
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Therefore, when collecting JSSl), (l5lT)D . (I5T3]) . KTEh . KTB and K2E\i . we 
deduce finally the desired bound, 

N 

Y.if^ . Hfj) 

- , (5.26) 

uniformly with respect to N and the orthonormal family {fj}- 

5.3. Upper bound. 

Coherent states for the curved boundary. 

Let a G (identified with dQ through a parametrization of the bound- 

ary). Let (/> = be the gauge function from Lemma [231 With $ = 
being the coordinate change near the boundary, let = B{^~^{a,0)) i.e. 
the magnetic field at the boundary point parameterized by a. Let finally, 
X G C°°(M) be a positive, smooth, have suppx C M+, and be normalized 
such that /igX^(s) ds = 1. Define Xh{s) for s G [0, [ to be 

and extend Xh to be a |5r2|-periodic function. Setting 
fj{{s,t);h,a,^) := 

(x)'^'""^'"^^^' (V?''^) ^'''^'^'^'^('-''^Mt), (5.27) 
we get — by the coordinate transformation <1> — the following function in ^l{tQ), 

fjix;h,a,0=fjiHxy,h,a,0. (5.28) 

We will consider fj as a function on all of O by extension by zero to Q\^}{tQ). 
Let n^"'^(/i, cj, .^) be the operator with integral kernel 

n^^°^(/i,f7,e)(x,x') = fj{x;h,a,Ofjix';h,a,0- (5.29) 

In terms of the projectors constructed in (|4.9p . n^°'^(/i, o", .^) is expressed as 
follows, 

n^^^'i(/i,a,0 = ^Pi,hU^'e-"^-/^Xhi- - cT)n, (^, B^;^)e"^-/^Xh{- - ct)C/* Vi,h • 
Then 

tr[n^^'^'i(/i,a,0] = / \fj{x;h,a,0\^dx 
Jn 

= ff \fj{{s,t);a,()\\l-tkis))dsdt 

< {1 + T{h)\\k\U II ^\u^{^t;i)\\{h)-W^)dsdt 

= (l + r(/i)||A;||oo). (5.30) 



m(e) 



B{x)_ 
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Also, using Lemma [231 

<(l + e) / \{hV - iB^Ao)e''^fjih,a,0\'^ dsdt 
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+ Ce"^ (r(/i)4 + /i2) J \f^{h,a,0\^dsdt 



+ 



T{hy 

/i2 



dsdt 



< (1 + e)iii{i)hB, + C — + e-^{T{hf + h 



(5.31) 



T{hf 

Let M{h,a,^,,j) be a function with < M < 1, and write 

, = g//M(...,e„)nr('...0^^. (5.32) 

Then, clearly < 7 as an operator on L'^{Q). We will prove that 

7 < (1 + l|A:||oor(/i)) . (5.33) 
Consider g G L'^{n) and define gg = gl{dist{x,dQ)<T{h)}- Then 

ialig) = {goliga), 

so we may assume that g = gg, i.e. that suppr^ C {x : dist(3;,50) < t(/i)}. 
We calculate, 



5(S) t)fj{s, t; h, a, ^)(1 — tk{s)) dsdt 



dadS, 



(5.34) 



We estimate from above by replacing JjM x | • p by // 1 x | • p in the above 
expression. Then we use the fact that Uj is an orthonormal basis of 
L^(]R+) for all ^. Next we evaluate the ^ (Fourier) integral, which becomes 



i^{s-s')./B:jh ^ 2-K^/h/B~„5{s - s') 

After inserting these two results, (I5.34p becomes 

{9\l9)<! I I \gis,t)\\l-tkis))^ 
Ja£dn Jsedn Jtm+ 

s - 
~{h) 



(5.35) 



T{h)-^X^{^-TSr)'^ih{'t)dsdtda. (5.36) 
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We do the cr-integration first. The normahzation of x imphes that the result 



IS 



{g\l9)< / \g{s,t)\'{l-tk{s)y^l^l,,it)dsdt 

< (l+r(/i)||A:||oo) / \g{x)\^dx. 
Jn 



That finishes the proof of (|5.33|) . 

Let K > 0. We choose the function M = ■ l{j=i}, with being 
the characteristic function of the set 

{(^, G X M : A _ ^,,(^) > 0, 1^1 < K} . 

Let 7^ be the corresponding density matrix. We calculate, 

M^(/i, a, h, a, 0) - bh\\M-, h, a, Ohnn)) 

< JJ M''{h,a,0[{l + e)fii{0hB^-bh{l- 



+ C 



/-l^^l dad£ 
-K Jo 



\B\ 



h{e\\B\\Lo.(^an) + r{h)\\k\\oo)m2K^^^-j= 



L°°{dn) 



C 



+ e-\T{h)^ + h^) \dn\2K 



\B\ 



The choice 

yields an error term which is controlled by CKh^^^. Therefore, 

-~^J kI B'^'iw^-^'M+dad^ + CKhy^ (5.37) 

Since K can be chosen arbitrarily large, this implies the upper bound in 
Theorem 11.11 

6. Proof of Theorem 11.41 

6.1. Lower bound. We follow the strategy of the proof of Theorem 11.11 
Let H = Ph^Q — bh — ah?/"^ . We seek a uniform lower bound of 
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for any orthonormal family {/i, /2, • • • , /n} in D{H). Using the notation of 
Section [5l we write (see (|5.6p ). 



N 2 N 

Y,{fj , Hfj) = J2 T^^fj > i^k,h{H - Vh)i'k,hfj) . (6.1) 

j=l k=l j=l 

The term corresponding to = 1 is a boundary term and is estimated exactly 
as in Section [5] (see also (|5.26p ). the extra term a/i^/^ being included in the 
error term. Notice that the localization error /i^r~^ = h^^^ ^ h^^'^ , so the 
a/i^/^-term can be incorporated without change to the argument. Therefore, 
the following lower bound holds, 

N 



27r 



/ b^/^ _d^ds{x)-0{h^/^^). (6.2) 

JdUxR 



Next, we estimate the interior term corresponding to = 2. Let us in- 
troduce a localization function x ^ ^"^(5(0, 1)) such that ||x||l2(i[52) = 1. 
Define 



1 / 



X — z 



2 



where C{h) g]0, 1[ is to be chosen later. Let us notice that 

Xhix;z)=0 VxGf], VzGJ] + 5(0,C(/i)). 
Then we get, using the standard localization formula, 

N 

Y.{fj,i^2MH-Vh)iP2,hfj) 

N 

^{fj ) i^2,hXhix; z){H - Wh)ip2,hXh{x; z)fj) dz , (6.3) 

J=l 

where Wh{x) = Vh{x) + /i^ f^^ l^^xXix] z)\^dz. 

The term corresponding to Wh will be estimated using a Lieb-Thirring 
inequality. Let 6 g]0, 1[ to be chosen later, and write, 

{fj , ip2,hXh{x; z){H - Wh)'ip2,hXh{x; z)fj) 

= (1 - ^){fj 1 ^{x; z)H(p{x; z)fj) 

+ 6{fj , (^(x; z){H,^^2 - 6-'WhMx; z)f,) , (6.4) 

where 

Whix) = ln{x) + T4(x) + j^^ |V,x(x; z)\^dz^ , 

iph{x;z) = 'il)i^h{x)xh{x; z) , 
and i^6R2 = — (/iV — i6Ao)^ — 6/1 is the self-adjoint operator acting in L^(M^). 
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As explained in (15.51) . the variational principle yields (i.e. Lemma [2?2l) . 

N 

To estimate the last term, we use the Lieb-Thirring inequality of Lemma 12.1 
Thus we get, 

TV 

^ ^{fj ' z){Hh,R2 - 6-^Wf;^)ip{x] z)fj) 



> 



r{h) V r{hY) Q{hY V 

Recall that we have already fixed a choice of 

T{h) = /l^/S . 

So, for the other parameters, choosing 

we get that the term on the r.h.s. above is of order h^^^ = o(/i^/^). Integrat- 
ing with respect to z € i7 + -6(0, C,{h)), we get, 

S / '^{fj , ^{x; z){Hh^^2 - 5 ^Wh)(p{x; z)fj) dz 
'' i=i 

>-C\9.\h^'^ = o{h^/^). (6.5) 

Now we calculate the leading order term. Recall the family of projectors 
(np)p>i on the Landau levels introduced in Subsection 14.11 We write, 

N 

^{fj > ^{x; z)Hip{x; z)fj) 
oo N 

p=i j=i 

oo N 

= E E^^J- ' ^^(^5 - m - ah^/^)U^{h, bMx; z)f,) 

p=i j=i 

N 

> -[a] + /i3/2 ^(/. , <^(x; z)Uf{h, bMx; z)fj) . 
The operator ip{x; z)Ili {h,b)ip{x; z) is positive. Therefore, 

N 

Y^ifj , ^{x;z)Uf{h,b)ip{x;z)fj) < tr {ip{x; z)Uf {h,b)ip{x; z)) . 
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Using (lO) . 

tr ((^(s z)Uj:{h, bM-; z)) = ^ vH^; z)dx<l. 

Hence, upon integrating w.r.t. z G + B{0, C{h)), we conclude the following 
lower bound, 



il-5)Y,{fj,V>ix;z)Hip{x-z)f,)> ^(1-5) / / ^\x;z)dxdz 



> -^^^[aU - Oih^/^) . (6.6) 



Combining (I6.1|1 - (I6.4|) . (|6.5p and (|6.6p . we get the following lower bound, 
which is what we desire to prove. 

6.2. Upper bound. We just construct a trial density 7 = 7int + 7bnd and 
estimate tr(i/7). We take 7int to be 

and 7bnd exactly as given in (I5.32j) . 

Notice that 7bnd and 7int act as direct sums since their integral kernels 
have disjoint support. 

By calculating tr (^^7) = tr (Hji^t) + tr (//7bnd) we will get the desired 
upper bound. The calculation of tr(ff7bnd) has already been carried out in 
Section [531 In order to calculate tr (Hji^t), we define, 

ip{x; z) := ^2 ( ^^^^TTW-^ ) Xh{x; z) , 7(2) = ip{-, z)U^{h, b)ip{-, z) . 



4/i3/8 



Then, 



tr(i77int) = j tT{H^{z))dz. 

Since II^^ is a projector, it follows that, 

tr(i?7(z)) = ii [Hip{x]z)Ii[{h,h)ip{x;z)) 

= tr (l[{{h,h)ip{x]z)Hip{x]z)Ii{{h,b)) . 
Using the following localization formula: 

{Hxf , Xf) = Re ixHf , xf) + h'W |Vx| V f , 

yx^c^m, y f€D{H)^ 

we deduce that 

tr {Uf{h,b)ip{x;z)Hip{x;z)Uf{h,b)) = 

Re[tr {uf{h,b)ip{x;zfHUj:{h,b))] + tr {uf{h,b)Vh{x;z)n^{h,b)) 



(6.7) 
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where 

Vh{x-z)=h'\Vip{x-z)\^ . 
Using that the trace is cychc, we get 

tr {Yl\{h,b)^{x]zfHIi\{h,b)) =tr {^{x; z)HYl\{h,h)^{x; z)) , (6.8) 

(6.9) 



tr {li{{h,h)Vh{x-z)Il{{h,h)) =tr [^VhIi{{KhUVh) . 

Now HUf{h,b) = -ah^/'^Uf{h,b). Therefore, we get, 

tr {ip{x; z)Uf{h, b)HUf{h, b)ip{x; z)) 

= -a/i^/^tr {ip{x;z)U^{h,b)ip{x;z)) . 

Notice now that, (recall that ip{x; z)Ili{h, b)ip{x; z) is a kernel operator), 
tr ((^(x; z)nf'(/i, 6)99(3;; z)) dz 



2'Kh 
b 

27rh 
b 



-02 



dist(x, d^) 



4/i3/8 



dx 



101 



/ dist(x,Oil)^ ^ 
V 4/i3/8 

, / dist(x, dfi) 

4/.3/8 



dx 



dx 



The function ipi (^^^^^^^^ being supported in $7(4/1^/^), its integral 

becomes small of the order 0{hP^^). Therefore, coming back to (|6.8p . we get 
finally. 



tr {ip{x;z)U^{h,b)HU^{h,b)ip{x;z)) dz 

We need next to estimate the trace (|6.9p . Actually, 
tr (^,U^{h,b)y^Vh^ 

Vh{x) dx 



ab\i}\ 

27T 



/,l/2 + 0(/i5/8). 



(6.10) 



2irh 

<^/ Iivur + 



dx 



Choosing ({h) = /i^/i^ then coming back to (|6.9p . we get, 
tr {ui{h,b)Vh{x;z)Uf{h,b))=0{h'>/^). 



(6.11) 
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We integrate ([O) w.r.t. z G + B{0,C{h)) and we substitute (fOOl) and 
(|6.11|) in the resulting formula to get, 
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